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Nomenclature 


In this paper, an integrated aerodynamic/dynamic 
optimization procedure is used to minimize blade weight and 4 
per rev vertical hub shear for a rotor blade in forward flight. 
The coupling of aerodynamics and dynamics is accomplished 
through the inclusion of airloads which vary with the design 
variables during the optimization process. Both single and 
multiple objective functions are used in the optimization 
formulation. The ’Global Criteria Approach’ is used to 
formulate the multiple objective optimization and results are 
compared with those obtained by using single objective 
function formulations. Constraints are imposed on natural 
frequencies, autorotational inertia and centrifugal stress. The 
program CAMRAD is used for the blade aerodynamic and 
dynamic analyses and the program CONMIN is used for the 
optimization. Since the spanwise and the azimuthal variations 
of loading are responsible for most rotor vibration and noise, 
the vertical airload distributions on the blade, before and after 
optimization, are compared. The total power required by the 
rotor to produce the same amount of thrust for a given area are 
also calculated before and after optimization. Results of this 
study indicate that integrated optimization can significantly 
reduce the blade weight and the hub shear as well as the 
amplitude of the vertical airload distributions on the blade and 
the total power required by the rotor. 
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chord 

natural frequencies of first four coupled elastic 
modes 

constraint function 

principle radius of gyration 

radius of gyration about reference axes 

area solidity 

th 

nonstructural weight per unit length at j node 
reference axes 

- i_ .th 

distance from the root to the center of the i 
segment 

box beam cross sectional area 
autorotational inertia 
power coefficient 

Young’s modulus 
bending stiffnesses 

objective functions 
factor of safety 
blade vertical airload 

amplitude of 4/rev vertical hub shear 

torsional stiffness 

principal area moments of inertia about reference 
axes 

length of l segment 
,.th 

mass of i segment 

number of nodes 
number of constraints 
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NDV 

number of design variables 

NSEG 

number of blade segments 

R 

blade radius 

W 

blade weight 

w. 

. . . .th 

weight of i segment 

W NS 

nonstructural weight 

W NS. 

nonstructural weight of i th segment 

w s ' 

structural weight 

w s. 

1 

structural weight of i* 1 segment 

Y 

nondimensionalized radial distance 

a 

prescribed autorotational inertia 

X 

taper ratio 

♦i 

i ^ design variable 

p i 

weight density of the 1 th segment 


advance ratio 

0 . 

i 

centrifugal stress in i* 1 segment 

o 

max 

maximum allowable stress 

V 

blade azimuth angle 

n 

rotor speed 


Subscripts 

r root value 

t tip value 

L lower bound 

U upper bound 

Introduction 

In recent years structural optimization has become a 
practical tool which can expedite mechanical design. An 
extensive amount of work has been done in developing design 
optimization procedures to bring the state of the an to a high 

level 1 ' 5 . While these techniques have received wide attention 

for fixed wing aircraft 1 , they are less well known in the rotary 
4-5 

wing industry . In the past the conventional blade design 
process was controlled mainly by the designer's experience 
and by trial and error. Today with improved understanding of 
helicopter analyses and efficient optimization schemes, 
attempts are being made to apply design optimization 
techniques and eliminate the expensive man- in-the -loop 
iterations. 

The helicopter rotor blade design process is 
multidisciplinary in nature and involves a merging of several 


disciplines, including structures, aerodynamics, dynamics and 
acoustics. For example, the blade design must satisfy 
specified strength criteria and be damage tolerant. The blade 

must be aeroelastically stable^ and the aerodynamic 

g 

performance should be optimized . The noise levels generated 
by the rotor which arc a function of local Mach number and 
9 

airloads must be minimized . 

An important criterion in rotor blade design has been 
reduced vibration. For a helicopter in forward flight, the 
nonuniform flow passing through the rotor causes oscillating 
airloads on the rotor blades which arc translated into 
undesirable vibratory shear forces and bending moments at the 
hub. Vibration alleviation therefore plays a major role in the 
10 

rotor blade design . In the conventional design process this 
has been accomplished by postdesign addition of 
lumped/tuning masses which may cause significant increases 
in the blade weight. Hence an important design consideration 
is to obtain an optimum blade design for reduced vibration 

without incorporating a weight penalty^* 11 

In most of the woik 8 ’ 1115 where optimization techniques 
were used to address blade design, attempts were made to 
satisfy certain design requirements and criteria related to a 
single discipline. The blade design was treated as a series of 
nearly independent design tasks with little consideration of the 
couplings and interactions between disciplines. For example, 
only blade aerodynamic requirements were considered in Ref. 
8. In Refs. 1 1 and 12 only dynamic design requirements were 
considered in optimum rotor blade designs. Taylor addressed 
the problem of vibration reduction through modal shaping in 
Ref. 11. In Ref. 12 Peters’ approach was to place the natural 
frequencies away from the critical frequencies and thereby 
reduce the hub shear. Blade structural requirements only were 
considered in Ref. 13. In Refs. 14 and 15 minimum weight 
designs were presented for articulated rotor blades with 
rectangular and tapered planforms but in the absence of 
airloads. 

The necessity of merging appropriate disciplines to 
obtain an integrated design procedure is being recognized. 
With improved understanding of helicopter analyses, it is now 
possible to include the couplings between the disciplines in the 
optimization procedure. Some initial investigations at partially 
integrating some of these disciplines arc presented in Refs. 6 
and 16-19. In Ref. 6 the dynamic and aeroelastic requirements 
were integrated and in Refs. 16-18 the dynamic design 
requirements were coupled with prescribed airloads in the 
analysis. Ref. 19 presented a formulation of the integrated 
procedure involving dynamics and aerodynamics. 
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Currently at NASA Langley Research Center, there is an 
effort to integrate various disciplines in the rotor blade design 
20 

process . The present paper is a pan of this effort. The 
purpose of the present work is to take an initial step in 
integrated aerodynamic/dynamic optimum design of a rotor 
blade by coupling two important technical disciplines - blade 
aerodynamics and dynamics. Design variables affecting both 
aerodynamics and dynamics are used and coupled 
aerodynamic and dynamic blade analyses are involved. The 
coupling is incorporated by including blade airload 
calculations in the analysis. The airloads change during the 
optimization process as a function of the design variable 
values. A trim analysis is pan of the procedure and the blade 
is trimmed at each intermediate design step of the optimization 
process. This paper concentrates on the low vibration and the 
low blade weight aspects of the design. The objective of the 
work is to reduce the. blade weight and the 4 per rev (4/rev) 
vertical hub shear, which is the critical shear force transmitted 
to the hub by a four-bladed rotor. When both the weight and 
the shear are simultaneously minimized a multiple objective 
function technique is needed. A method known as the ’Global 
21 

Criteria Approach’ is used in the present work. The 
optimization procedure developed is applied to a blade test 
problem. Results obtained with a multiple objective function 
formulation are compared with two single objective function 
formulations - one with the blade weight as the objective 
function and the other with the 4/rev vertical shear as the 
objective function. 


Problem Description 

The design goal is (0 minimize the weight and the 4/rev 
vertical hub shear of a rotor blade in forward flight. The 
constraints are ’windows’ on the coupled flap-lag natural 
frequencies to prevent them from falling into the critical 
ranges. These windows are selected by careful consideration 
of the blade frequency versus rotor speed diagram and the 
associated frequency coalescence points. Constraints are 
imposed on frequencies of the first four coupled elastic modes. 
A prescribed lower bound on the blade autorotational inertia 
and a maximum allowable upper bound on the blade 
centrifugal stress are also included as constraints. 


is the torsional stiffness, k is the radius of gyration and w. is 

the nonstnictural weight at the j th node. The subscript V 
refers to values at the root and N denotes the number of nodes. 
The lumped nonstructural weight at node 1 is not used as a 
design variable since it involves contribution from the hub. 

Optimization Formulation 

Three optimization formulations are studied. The first 

formulation minimizes the blade weight W, the second 

minimizes the 4/rev vertical shear F and the third minimizes 

z 

the weight and the 4/rev shear simultaneously. To minimize 
the blade weight and the 4/rev vertical shear simultaneously, 

21 

the Global Criteria Approach is used. Using this method the 
optimum solution is obtained by minimizing a ’global 
criterion’ F(0) where 


F(0) 


/ W( $) -W(**)\ 2 /F z ( +) - F z ( ♦*) \ 2 

V W(#*) / \ F 2 ( ♦*) / 


subject to 


( 1 ) 


gj($) < 0 j » 1, 2, .... NCON. 


The global criterion is thus defined as the sum of the squares 
of the relative deviations of the individual objective functions 
W and F^ from their respective individual optimum values. 

♦ 

The design variable vector 0^ is obtained by minimizing the 

single objective function W(0) subject to the set of constraints 

* 

g(0) and the design variable vector 0^ is obtained by 

minimizing F^ subject to the same set of constraints. In Ref. 

19 this approach gave good results when minimizing weight 
and centrifugal stress for a blade in vacuum. However, in the 
present case, due to the highly nonlinear nature of the above 
function the following objective function is used 


A linear taper is allowed along the blade planform (Fig. 
1). The blade taper ratio, X, is expressed as follows : 


X = c / c 
r t 


where is the root chord and is the tip chord. The design 

variables are the following: El , El , GJ . X, c , k , and w. 

xx zz r r r j 
r r 

( j=2 N) where El and El are the bending stiffnesses, GJ 


xx 


zz 


F(0) WF(0) 

where F(0) is given by Eqn. 1. 

The three formulations are summarized in Table 1. In 
th 

this table represents the k natural frequency where k=3,..,6 

(the first two modes are assumed rigid body modes). The 
subscripts L and U denote lower and upper frequency bounds, 
respectively. The quantity a is the prescribed lower bound on 
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(5) 


the autorotational inertia AI, o. is the centrifugal stress on the 


i blade segment, a is the maximum allowable blade stress, 
max 

FS is a factor of safety and NSEG denotes the number of blade 
segments. 


I 

xx 


k 

xx 


I 

zz 



( 6 ) 


Blade Analytical Reprraptalion 

In this section, expressions for cross sectional properties 
are presented for the blade model (Fig.l) used in the 
optimization process. Since the detailed design geometry of 
the blade cross section is not pan of the current work, the 
cross sectional area and the weight of the blade are formulated 
in terms of the design variables used. The total blade weight, 
W, has two components - the structural component W ^ and the 

nonstructural component W^. The blade weight is: 

NSEG 

W- £ (W S .W NS , (2) 

i =1 i i 

where and are the structural and nonstructural 

i i 

weights of the I th segment respectively. The structural weight 
of the i * blade segment, , is 


where A denotes the cross sectional area and k and k are the 

xx zz 

radii of gyration about the z and x axes respectively. The 
expression for A, obtained by summing Eqns. 5 and 6, is given 
by 

A * (EI^+Ey/Ek 2 (7) 


where k (principal radius of gyration) is 


,2 2 
k + k . 
xx zz 


( 8 ) 


Assuming that the blade stiffnesses are contributed by the 
blade structural component only and recalling that the inertia I 
4 

* L , the distribution of the stiffness EI^ along the blade radius 
is expressed in terms of its root value as follows 



i 


A. p.L. 

l l t 


(3) where R is the blade radius* y is the radial location and c(y) is 
the chord distribution for a linear taper given by 


where A., p. and L. denote the segment cross sectional area, 

lit 

the weight density and length respectively. Since the 
nonstructural weights are specified at the nodes in the test 
problem the segment nonstructural weight is calculated by 
averaging the values of these weights at the two associated 
nodes as follows 


W 


NS. 


L i‘VVi > /2 


<4) 


where w. denotes the nonstructural weight per unit length at 


c(y) 




(10) 


Similar expressions are derived for El^Cy) and GJ(y). 


The autorotational inertia (AI) of the blade is calculated 
from the segment weights as follows 


NSEG 

A,. £ w,V 

t *1 


( 11 ) 


the i lh node (Fig. 2). 

The expression for the blade cross sectional area 
contributing to the blade structural weight is developed in 
terms of the radii of gyration and the stiffnesses as follows. 
Assuming that the blade is made of isotropic material the 
moments of inertia of the blade structural component about the 
z and x axes (Fig.l ) are given by 


where W. denotes the total weight of the i^ segment and y. 

denotes the distance from the root to the center of the 

segment. The centrifugal stress in the i^ segment is calculated 
as follows 
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o. 

1 



A . 

i 


(12) 


where o. is the stress due to centrifugal forces in the i 
segment and Q is the rotor speed. 


A na l ys is 


The optimization procedure with integrated 
aerodynamics and dynamics is described next. The procedure 
uses Comprehensive Analytical Model of £otorcraft 

22 

Aerodynamics and dynamics (CAMRAD) for the blade 
aerodynamic and dynamic analyses. The basic algorithm used 
for optimization is the method of feasible directions 

implemented in the program CONMIN 23 along with an 
approximate analysis. 

Agr^ ynamitAnalKsia 

CAMRAD is used to calculate the section loading from 
the airfoil two-dimensional aerodynamic characteristics. It 
uses the lifting line or blade clement approach and has 

24 

corrections for yawed and three-dimensional flow effects 
The program also has the provision for including unsteady 
aerodynamics. It offers two broad categories of trimming - the 
free flight case and the wind tunnel case. The wind tunnel 
trim option is selected for this work since the model used in 
this study is a wind tunnel model of a rotor. The trim option 
consists of trimming the rotor lift, drag and the flapping angle 
with collective pitch, cyclic pitch and shaft angle. It is, in fact, 
a coupled aerodynamic and dynamic procedure. The blade is 
trimmed for the specified flight condition at each step of the 
design process. The procedure incorporates the change in 
planform caused by design variables during optimization. 


Dynamic Analysis 

The modal analysis of the blade is also performed in 
CAMRAD which uses a modified Galerkin approach to 
calculate the mode shapes and frequencies. The method works 
better than the Rayleigh Ritz approach in the presence of large 

24 

radial variations of the bending stiffness . Using the airloads 
information obtained from the aerodynamic analysis the 4/rev 
vertical shear is calculated. 


Optimization Implementation 

The optimization procedure is shown in Fig. 3. The 
process is initiated by identifying the blade preassigned 


parameters which are held fixed during optimization. The 
next step is to initialize the design variables and perform the 
blade structural analysis to calculate the blade properties, the 
centrifugal stress and the autorotational inertia. The 
aerodynamic and dynamic analyses are performed next using 
CAMRAD as described in the previous section. A sensitivity 
analysis is part of the procedure and consists of evaluations of 
the derivatives of the objective function and the constraints 
with respect to the design variables. Analytical derivatives are 
used for the weight, autorotational inertia and centrifugal 
stress. Forward finite differences are used for the derivatives 
of the hub shear and the frequencies. Once the sensitivity 
analysis is completed optimization is performed. 

The optimization process requires many evaluations of 
the objective functions and constraints before an optimum 
design is obtained. Therefore an approximate analysis is used 
along with CONMIN to save computational effort. In the 
present work, the objective function and constraints are 
approximated using a piecewise linear analysis that consists of 
linear Taylor series expansions for the objective function and 
the constraints. Using CONMIN, along with the approximate 
analyses, updated design variable values are obtained. The 
process continues until convergence is achieved. Convergence 
is checked on the objective function value over three 

consecutive cycles using a tolerance of 0.5 X 10 

TsslProbtem 

The baseline blade model selected is a wind tunnel blade 

25 

model of the Growth Black Hawk (GBH) which is a four- 
bladed rotor. The blade is articulated, has a rigid hub, a linear 
twist distribution and a distribution of three airfoils along the 
span. The GBH blade model has rectangular planform for 80 
percent of the blade radius and a taper ratio of 3 for the 
remainder of the blade. The initial CAMRAD model of this 
blade used 14 aerodynamic segments and 38 structural nodes. 
Some changes have been made to this model to make it more 
suitable for the optimization study, e.g., the structural 
discretization is modified to reduce the total number of design 
variables. This modified blade model will be referred to as the 
Modified Growth Black Hawk (MGBH) and has the following 
modifications 

1) the blade is tapered linearly from root to tip, 

2) eight structural nodes are used, 

3) the same airfoil is used for all aerodynamic segments. 

These modifications produce only minimal changes in the 
aerodynamic and dynamic performance behavior of the blade 
from that of the GBH blade. This has been verified by 
performing integrated aerodynamic/dynamic analysis of the 
two models using CAMRAD. Table 2 contains the parameters 
of the MGBH and the GBH blades which are fixed during 
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design optimization. It should be noted that the GBH model 
rotor has achieved mission requirements in excess of the 
requirements used in the optimization problem. Therefore the 
MGBH rotor blade also would satisfy the requirements to 
begin with though not efficiently. The MGBH blade is used as 
the starting blade design for the optimization study in the 
present paper* with a taper ratio of 1 , and will also be called 
the ’reference blade.* 

R«ultS-and-Disc»s§iQn§ 

This section of the paper presents results from the three 
formulations summarized in Table 1. The aerodynamic 
analysis of the blade is performed using CAMRAD with 
uniform inflow. The blade is in forward flight with an 
advance ratio, \i * 0.30. The frequency bounds summarized in 

Table 3. A minimum value a of 19.748 lb-ft^ is used for the 

autorotational inertia constraint. An allowable stress o of 

max 

2.0 x 10^1b/ft^ and a factor of safety FS = 2 are used. 

Results of the three formulations are presented. The 
effect of optimization on the blade dynamics is assessed by 
studying the nature of the first four elastic modes associated 
with the frequencies, - f^, constrained during optimization. 

The effect of optimization on airload distributions is discussed 
and the total power required by the rotor investigated. In all 
the optimization results obtained, the center of mass of the 
blade is forward of the quarter chord point. This is due to the 
variations of the stiffnesses assumed along the blade radius 
which do not allow large perturbations of the center of mass. 
Table 4 summarizes the results obtained using the three 
formulations. Column 1 presents the reference blade data and 
columns 2*4 presents the corresponding information for the 
optimum blades. In the formulations presented here optimum 
results have been obtained within seven to ten cycles. 

Weight Minimization (FormulationJi 

The design obtained using blade weight as the single 
objective function is presented in column 2 of Table 4. As 
shown in the table, it is possible to reduce the weight of the 
blade by 7.5 percent from the corresponding reference blade 
value while satisfying all the imposed constraints. After 
optimization, the fourth natural frequency, f^, is at its 

prescribed upper bound and the autorotational inertia 
constraint is at its prescribed lower bound. The optimum blade 
has a taper ratio of 1.17 and the values of the stiffnesses, El 

r 

and El , are increased from their corresponding reference 
zz 

r 

blade values. However, the value of the torsional stiffness at 
the blade root, GJ^ remains unchanged after optimization 


since no constraints are imposed on the blade torsional 
frequencies. The values of the design variables and c are 

also reduced after optimization. The reduction in the value of 
c^ reduces the blade area solidity as shown in the table. The 

value of the 4/rev vertical shear reduces by 40.6 percent from 
the reference blade value. 

Hub Shwr Minimization (for m ulation p 

The design with 4/rev vertical hub shear as the single 
objective function is summarized in column 3 of Table 4. The 
hub shear is reduced by 85.6 percent compared to the 
reference blade value. The third coupled natural frequency, , 

reaches the prescribed lower bound and the autorotational 
inertia constraint is at its prescribed lower bound. The 
optimum blade is tapered with the taper ratio 1.18. The values 
of El and El are higher than their corresponding reference 

r r 

blade values. The value of GJ^ remains unchanged during 

optimization for the reason described previously. The values 
of k^ and c^ are once again reduced after optimization and so is 

the blade solidity. The total Made weight is reduced by 8.5 
percent. The reduction is more significant in this formulation 
than in the previous formulation where the blade weight is the 
objective function. This can be explained as follows. The 
reduction of hub shear reduces the amplitudes of the blade 
radial airload distribution, as will be shown later, which 
reduces the blade stiffnesses and thereby reduces the blade 
weight. 

Combined — Weight — and Hub Shear Minimization 

(Formulation 3) 

Results of the simultaneous minimization of the blade 
weight and hub shear using the Global Criteria Approach are 
summarized in column 4 of Table 4. The 4/rev hub shear is 
reduced by 77.6 percent and the blade weight by 10.6 percent. 
The constraints are all satisfied and the autorotational inertia 
constraint is greater than its prescribed lower bound. The 
optimum blade has a taper ratio of 1.33, which is larger than in 
the other formulations. The optimum values of the bending 
stiffnesses El and El are higher than their respective 
r r 

reference blade values and the value of the torsional stiffness, 



therefore in the value of blade solidity are the largest. 


Additional Comparisons of Designs 

The distributions of the nonstructural weight/unit length 
for all three formulations, before and after optimization, are 
shown in Fig. 4. In all formulations the nonstructural weight 
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increases outboard after optimization due to the presence of 
the autorotational inertia constraint. The trend is more 
significant in formulation 2 (hub shear) than in formulation 1 
(weight) and is most significant in formulation 3 (weight and 
hub shear). It is interesting to note that in the hub shear and 
the multiple objective function formulations, the maximum 
values of nonstructural weights occur at the next to last node 
and not at the last node as in the weight optimization 
formulation. This could be due to some changes in the 
coupled mode shapes effected by modified stiffness 
distributions when the hub shear is minimized either 
separately or simultaneously with the weight. Fig. 5 presents 
comparisons of optimum weight and vertical shear from all 
three formulations. Fig. 5a compares the blade weight and 
Fig. 5b compares the 4/rev vertical shear. As shown in Fig. 5 
and Table 4 the Global Criteria Approach provides the lightest 
blade structure with a hub shear reduction that falls between 
those obtained from the single objective function 

optimizations. This is contrary to the intuitive belief that the 
use of a multiple objective formulation should yield solutions 
lying between those of the single objective formulations. In 
other words, the blade weight obtained by simultaneously 
minimizing weight and hub shear might be expected to be 
higher than that obtained from weight minimization and the 
hub shear obtained should be higher than that obtained from 
hub shear minimization. However, this is only true if the 
objective functions are monotonically increasing functions of 
the design variables. This is not true in the present case since 
for example, the blade weight can decrease with an increase in 
taper ratio and the hub shear is a very complicated and 
nonmonotonic function of the design variables. 

Effect of Optimization on Mode Shapes 

There was a concern that the nature of the mode shapes, 
corresponding to the frequencies constrained in the 
optimization procedure, might change as a result of 
optimization. Therefore the mode shapes for the first four 
elastic modes are plotted before and after optimization in Figs. 
6-9. Pan (a) of each figure shows the flapping components 
and pan (b) shows the lead- lag components of the coupled 
modes. The figures indicate that there is a change in the 
magnitude of the modal deflections after optimization and in 
most cases, the node points of the modal deflections are 
moved slightly outboard from the reference blade. However, 
no major changes in the nature of the mode shapes occur. 

Effect of Optimization on^lad^Aerodvnainics 

Effect on Airloads - Reduction of hub shear is achieved when 
either of the following occurs: blade natural frequencies move 
away from the excitation frequencies; mode shapes become 
more nearly orthogonal to the forcing function or blade 
harmonic airloads are reduced. It was of interest in the present 


work to investigate the extent to which the optimization 
process alters the distribution of airloads. 

The influence of the design variables on the section 

aerodynamic load normal to the plane of the rotor, F , is 

P 

assessed by comparing its radial and azimuthal distributions 
before and after optimization. Fig. 10 presents the radial 

distributions of F for the advancing blade at 90 degree, 
P 1 

for both the reference and the optimum blades. Fig. 10 
indicates that optimization reduces the amplitude of the 
vertical airload throughout the span. The reduction is largest 
in formulation 3. Note that the high negative value of F^ near 

the tip is caused by the high negative twist of 16 degrees 
(Table 2 ) in which the tip of the blade has a smaller pitch 
angle than the root. Blades with negative twist may carry 
negative lift outboard and positive lift inboard on the 
advancing side. It is to be noted that the overall thrust and the 
propulsive forces are the same for the reference and the 
optimum blades. 

In Figs. 11-13 the azimuthal distributions of the vertical 

airload, F , are plotted for three radial locations. The radial 
P 

locations selected are near the blade root, at the point of thrust- 
26 

weighted equivalent chord and near the blade tip. Fig. 1 1 
presents the azimuthal distributions of F^ near the blade root 

(Y = 0.25, where Y denotes the nondimensionalized radial 
location). Fig. 12 presents the azimuthal distributions at Y = 
0.75 which is the point of occurrence of the thrust- weighted 
equivalent chord for all the three optimum blades. Fig. 13 
shows similar distributions near blade tip (Y = 0.9875). These 
figures indicate that optimization reduces the amplitudes of the 
vertical airloads at all three radial locations. Further, 
formulation 3 produces the maximum reductions. It can 
therefore be concluded from Figs. 11-13, that the design 
variables selected for the optimization formulation have 
significant influence on the amplitude of blade radial and 
azimuthal distributions of airloads. To complete the 
investigation the vertical airload distribution 

nondimensionalized with respect to solidity, F /s, is plotted to 

P 

study the effect of optimization on the vertical blade loading 
for a given area. Fig 14 shows that the optimum blades 
produce a larger vertical force for a given area, F /s, without 

p 

stalling, i.e., at any radial location the optimum blades use the 
area more effectively. The best results are obtained using 
formulation 3. 

Effect on Power Required - Finally, the total power required 
by the rotor to provide the same amount of thrust is checked 
before and after optimization. Figs. 10-13 show that 
optimization not only reduces the amplitude of the vertical 
airload distributions but it also reduces the peaks and valleys 
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making the distributions more harmonic both radially and 
azimuthally. This should indicate a reduction in the total 
power required by the rotor. Therefore, the power coefficient 
is calculated before and after optimization. The 
nondimensional power coefficient C^/s is shown in Fig. 15 for 

the reference and the optimum blades. The figure indicates 
that the total power required to produce the same amount of 
thrust for a given area reduces significantly after optimization 
and the reduction is most significant in formulation 3. Hence, 
optimization improves not only the net rotor power 
requirement which should reduce since the solidity of the 
blade is reduced but it also reduces the power per unit area. 

i 

Cflpclusipns 

In this paper rotor blade designs are generated by use of 
an integrated aerodynamic/dynamic optimization procedure. 
Blade weight and 4/rev vertical shear are minimized in the 
presence of airloads under forward flight conditions using both 
single and multiple objective functions formulations. 
Constraints are imposed on the first four coupled natural 
frequencies (elastic modes only), the blade autorotational 
inertia and the centrifugal stress. The program CAMRAD is 
used for the blade dynamic and aerodynamic analyses and the 
program CONMIN is used for optimization. The Global 
Criteria Approach is used for the multiple objective 
formulation and its results are compared with those obtained 
from single objective function formulations. The optimum 
designs are compared with those of a reference blade. The 
optimization scheme used has been very efficient and 
optimum results have been obtained in seven to ten cycles. 

The following observations are made from this study. 
The Global Criteria Approach for formulating the multiple 
objective function optimization is very effective. The 
approach yields a design in which the blade weight is lower 
than that obtained by using the blade weight as the single 
objective function due to the nonmonotonicity of the objective 
functions with respect to the design variables used. The hub 
shear is somewhat higher than that obtained by using the hub 
shear as the single objective function. The optimum blades are 
all tapered in contrast to the reference blade which is 
rectangular. In all three optimization formulations studied the 
center of mass of the blade remains forward of the quarter 
chord point. Inspection of the airload distributions for the 
initial and the optimum designs indicate that optimization 
significantly reduces the amplitude and alters the distributions 
favorably to help achieve the goal of low vibration. The 
optimization procedure changes the amplitude of the modal 
deflections and in most cases moves the node points of the 
associated modes slightly outboard but no significant changes 
in the nature of the mode shapes occur. Calculation of the 
total power required by the rotor to produce the same amount 


of thrust for a given area indicates a significant reduction after 
optimization. 
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Table 1 Summary of the optimization problem 


Formulation 

1 2 

3 

Objective function 

W F 

WandF 

Constraints 
Design variables 

z z 

■-'A, 

k=3 6 

vv 1 

1 - AI/a 

FS o./o -1 (i=l,...,NSEG) 

i max 

El , El , GJ , X, c , r , w. 

XX zz 

r r 

r r r j 

(j=2,...,N) 


Table 2 Growth Black Hawk and Modified Black Hawk blade 
parameters 


Blade radius 

4.685 ft 

Number of blades 

4 

Rotational speed 

639.5 

Flap hinge offset/radius 

0.0534 

Inplane hinge offset/radius 

0.0534 

Maximum twist (at blade tip) 

-16 degrees 

Advance ratio 

0.3 


Table 3 

Blade frequencies (per rev) and bounds (windows) 


Reference 

Prescribed bounds 



lower 

upper 

f 3 

3.068 

3.05 

3.50 

f 4 

6.763 

6.50 

6.90 

f 5 

9.283 

9.25 

9.5 

f 6 

12.632 

12.5 

12.75 


9 



Table 4 Optimization results 


Formulation 


Reference 

1 

2 

3 

El (lb-ft 2 ) 

XX 

r 

10277 

11305 

10306 

11818 

El (lb-ft 2 ) 
zz 

r 

354 

389.42 

384.59 

402.15 

GJ r (lb-ft 2 ) 

261 

261 

261 

261 

k (ft) 

0.268 

0.192 

0.175 

0.176 

X 

1.0 

1.17 

1.184 

1.33 

c (ft) 
r 

0.450 

0.403 

0.295 

0.281 

f^ (per rev) 

3.068 

3.194 

3.271 

3.282 

f 4 (per rev) 

6.763 

6.853 

6.87 

6.815 

f 5 (per rev) 

9.283 

9.415 

9.25 

9.487 

fg (per rev) 

12.632 

12.750 

12.749 

12.512 

AI (lb-ft 2 ) 

19.748 

19.748 

19.748 

19.955 

4/rev hub shear Ob) 

0.160 

0.095 

0.023 

0.036 

Percent reduction in 
4/rev hub shear 

— 

40.61 

85.6 

77.6 

Blade weight Ob) 

3.408 

3.152 

3.120 

3.048 

Percent reduction in 
blade weight 

— 

7.51 

8.5 

10.6 

Blade solidity 

0.122 

0.098 

0.071 

0.062 



Fig. 2 Discretized blade model. 



Fig. 3 Flow chart of optimization procedure. 



Fig. 1 Simplified rotor blade model with linear taper. 
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4.0 r 


Weight (lbs) 


3.408 

! I 3.152 iim 



Rafarenca Weight Hub shear Weight & 

hub shear 

a) Blade weight 



hub ahaar 


b) 4/rev vertical shear. 


c Reference 



b) Lead-lag deflections. 
Fig. 6 Concluded. 


o Reference 



-1 L — — 4 * * * ■ 

0 0.2 0.4 0.6 0.8 . 1 

Nondlmensional radius, Y 


Fig. 5 Comparison of optimum weight and hub shear. 


a) Flapping deflections. 




Fig. 6 Deflections of first elastic mode. 


^ig. 7 Deflections of second elastic mode. 
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„ r c Reference 

.3 
.2 
.1 

Deflection 

0 

-.1 

-.2 

0 0.2 0.4 0.6 0.8 1 

Nondimensional radius, Y 

a) Flapping deflections. 



b) Lead-lag deflections. 



Fig. 9 Concluded. 



o Reference 



b) Lead-lag deflections. 


Fig. 8 Deflections of third elastic mode. 


Fig. 10 Radial distributions of vertical airload, F 


V = 90 degrees; p = .3. 



40 
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, v 

o Weight 
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o Hub shear 

Vertical 

A Weight & 
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Azimuth angle, y (degree) 


360 


a) Flapping deflections. 


Fie 9 Deflections of fourth elastic mode. Fig. ^ Azimuthal distributions of vertical airload, F ; 

P 

Y = 0.25; p = .3. 
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c Weight 
o Hub shear 



Fig. 12 Azimuthal distributions of vertical airload, F ; 

P 


Y = 0.75 ; h = .3. 
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.005 
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.002 


.001 


.0040 
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.0039 
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Hub shear Weight & 
hub shear 
9.6%* 
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* Percent reduction from reference blade 


Fig. 15 Comparison of nondimensional power coeffecient, 


C p /s. 


15 r 



Fig. 13 Azimuthal distributions of vertical airload, F ; 

P 

Y = 0.9875; p = .3. 



Fig. 14 Radial distributions of nondimensional vertical 

airload, F / s; v = 90 degrees; 4 = .3. 
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